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symbol for projectile particle in the reaction 1 + 2 — > anything 
symbol for target particle in the reaction 1 + 2 — > anything 
symbol for center of momentum frame, where p : + P 2 = 0 
symbol for lab (target) frame, where P 2 = 0 
denotes a variable x in the cm frame 
denotes a variable x in the lab frame 

denotes a variable x in the lab frame or denotes an invariant variable x 

symbol for pion (either 7r°,7r + , or tt~) 

symbol for nucleon (either proton or neutron) 

symbol for neutron 

symbol for proton 

Note: not to be confused with magnitude of momentum. Context will make clear. 

pion mass, GeV 

proton mass, GeV 

mass of particle i, GeV 

particle kinetic energy, GeV 

particle total energy, GeV 

magnitude of 3-momentum, p = |p|, GeV 

Note: not to be confused with symbol for proton. Context will make clear. 

speed of cm frame 

relativistic 7 factor, 7 = . 1 

1 ’ ' 07 ^ 

Lorentz invariant differential cross section, mb/GeV 2 
spectral distribution cross section, mb/GeV 
total cross section, mb 

single particle distribution in the projectile-nucleon cm frame 

total energy in cm frame, Ei ^ + E 2 cm = E cm , GeV 
pion momentum in cm frame, GeV 
maximum pion momentum in cm frame, GeV 
maximum pion angle in lab frame 

magnitude of 3-momentum of particle 1 in lab frame for reaction 1 
| P 1 1 lab) GeV 

total energy of particle 1 in the lab frame for the reaction 1 + 2 — >■ 

\/\p\uah \ rn t GeV 
transverse momentum, GeV 
4-momentum vector, = (E, p), GeV 
square of 4-momentum vector, (p ^) 2 = E 2 — p 2 = m 2 , GeV 2 
4-momentum vector of particle i, GeV 


+ 2 -> 3 + 4, 
3 + 4, 



Abstract 


Ranft has provided parameterizations of Lorentz invariant differential cross 
sections for pion and nucleon production in pion-proton collisions that are 
compared to some recent data. The Ranft parameterizations are then numer- 
ically integrated to form spectral and total cross sections. These numerical 
integrations are further parameterized to provide formulas for spectral and to- 
tal cross sections suitable for use in radiation transport codes. The reactions 
that are analysed are for charged pions in the initial state and for both charged 
and neutral pions in the final state. 

1. Introduction 

The development of the HZETRN code (refs. 1-4) which predicts space radiation 
environments, has reached the stage where meson interactions need to be considered. 
The lightest meson is the pion and is therefore the most important one from a radiation 
point of view. The next lightest meson is the kaon, and its interactions will be dealt with 
in future work. For the present, only pion interactions will be considered. 

Previous work (refs. 5-8) evaluated pion production in proton-proton collisions. The 
7 r° decays electromagnetically with a lifetime of r = 8.4 x 10” 17 sec that corresponds to a 
decay length cr = 25.1 mn. This short decay length means, from a radiation transport 
point of view, that it can be considered to decay instantaneously. Thus, for the HZETRN 
code, the 7 r° does not propagate, and its interactions with the target constituents do not 
need to be considered. 

On the other hand, the charged pions, tt + and 7r~ , decay weakly and therefore have a 
much longer lifetime and longer decay length. The lifetime is r = 2.6 x lCT 8 sec that cor- 
responds to a decay length cr = 7.8 m. This long decay length means that charged pion 
propagation needs to be included in HZETRN because the pions will interact with target 
materials before they decay. From this point of view, it is very interesting that the Weak 
interaction has an extremely important role in space radiation transport applications. 

The present work involves the interaction of these charged pions with the protons 
in the target material. The final states considered will be pion production and proton 
production. HZETRN requires spectral distributions in the lab (target) frame as well 
as total cross sections. Parameterizations of these cross sections will be developed as 
will parameterizations of the Lorentz invariant differential cross sections. (These latter 
parameterizations will be useful when developing angular distributions for later use in a 
3-dimensional transport code.) 



Analysis of the relevant pion reactions has been extensively discussed in reference 9. 
All possible Strong exclusive reactions were itemized, along with thresholds and mass 
values of the set of inclusive particles for the minimum threshold reactions. These results 
(ref. 9) are used in the present work, especially in doing numerical integrations. 

Appendices A-D contain the following information: appendix A presents data that 
relates to Lorentz transformations; appendix B has useful kinetic results; appendix C 
presents data about maximum pion momentum in the cm frame, and appendix D contains 
data regarding the maximum pion angle in the lab frame. 

The threshold energies are derived in reference 9. Because of the method used to 
parameterize cross sections (refs. 5-9), the threshold does not appear automatically in 
the parameterizations. When putting cross sections into a code like HZETRN, it is 
important to include a code statement that will set the cross section to zero when the 
energy is below threshold. 

2. Ranft Parameterizations of Lorentz Invariant Differential 
Cross Section for Pion-Proton Collisions 

Ranft (ref. 10) provides parameterizations of the single particle distributions, dp ~ dp | 
which are related to the Lorentz invariant differential cross section according to 

dPa _ d 2 N ^ 

7 i ^inelastic 7 7 l-U 

dp\\dp ± dp\\dp ± 

The total inelastic cross section (ref. 10) rises approximately as A 2 / 3 . For pion-nucleus 
collisions, Ranft (ref. 10) uses 


inelastic = 28.5A°' 73 (2) 

where A = 1 for the pion-proton reactions considered in the present work. This value 
of u inelastic is used in both of the following subsections for pion and nucleon (proton or 
neutron) production from pion-proton collisions. 

We require the Lorentz invariant differential cross section d pl c j E = E ^ where E = 
is the pion total energy and p = | | is the magnitude of the pion 3-momentum. The 

3-momentum volume element is 

dp 3 = p 2 dpdQ = 2'Kp 2 dps\n0d0 (in spherical coordinates) (3) 

= 2itdp\\p±dp± = 7rdp||(ip 2 (in cylindrical coordinates) (4) 
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We want the relation between -Afs; and , d *7 . Thus 

ap* j hi dp\\dp±. 


d 3 a 

dp*/E 


2.1 Pion Production 

This subsection considers the reactions 

7T T p — » 7T + X (7) 

Ranft (ref. 10) has parameterized the single particle pion distributions as 

for p\\ < 0 (i.e. 9 > 90°) (8) 

for p\\ > 0 (i.e. 6 < 90°) (9) 

where p± and p\\ refer to the final state pion produced and where s = E 2 m and E 2 = 
p‘f\ +P± + m n = P 2 + m \- (Note that Ranft incorrectly writes p\\ > 0 for both of the above 
formulas.) The Ranft parameters for pion production are listed in table 1. 

Table 1. Parameters for Equations (8) and (9) or (10) and (11) 
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( 6 ) 


Reaction 

CL\ 

a,2 

a 3 

a 4 

as 5 

Q,q 

d'j 

as 

7T + + p — » 7T + + X 

0.22 

7 

0 

5.7 

0.22 

5 

0.115 

5.7 

7T + + p — > 7T° + X 
7r"*~ T p — y 7T~ T X 
7r~ T p — y 7t° T X 
7T~ + p — » 7T + + X 

0.14-0.14/Vs 

13 

0 

5.7 

0.14- 0.14/Vi 

9 

0 

5.7 

7T~ + p — > 7T~ + X 

0.22 

7 

0 

5.7 

0.22 

5 

0.115 

5.7 


Substituting the Ranft parameterizations for 


dpN 

dp\\dp ± 


gives 


d 3 a 

dp*/E 


.2 r 


^inelastic 


27 r 


. tig 2 

a 5 exp(-4— P||) + a 7 


exp(-a 8 p ± ) 
for p\\ < 0 (i.e., 6 > 90°) 


(10) 
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I'K 


&2 2 \ 

“1 exp(— 4 — p||) + a 3 


exp(— a 4 pjJ 

for p|| > 0 (i.e., 9 < 90°) 


( 11 ) 


where E = and p 3 , p» and p± are also all final state pion variables. 


These parameterizations are compared to data in figures 1-10. Clearly the Ranft 
parameterization provides a poor fit to the data. 


2.2 Nucleon Production 


This subsection considers the reactions 


7T T p — > N + X (12) 

where N is either a proton (p) or neutron (n). 

Ranft (ref. 10) has parameterized the single particle nucleon (proton or neutron) 
distributions as 


d 2 N 

dpi | dp i 


~F ^1 T ^=p\\ + p± [ ex P ( ) + a 5 exp(-a 6 p_ L ) 


for p\\ < 0 (i.e. 9 > 90°) 


(13) 


and 

d 2 N 

— — - — = 0 for p\\ > 0 (i.e. 0 < 90°) (14) 

dp\ | dp± 

where p± and p\\ refer to the final state nucleon produced. The Ranft parameters for 
nucleon production are listed in table 2. 


However, these parameterizations cannot possibly be right as the data (ref. 11) show a 
large cross section for 9 < 90° and a smaller, non- zero cross section for 9 > 90°. Thus, for 
the rest of this section and in figures 10 - 15, the Ranft parameterization is herein defined 
as equation (13). Thus, we define the Ranft parameterization for nucleon production as 


d?N 

dp\\dp± 


CL\ 


, , a ‘i “3 2 

1 + + 7 P » 


P . L 


exp(-a 4 p^) + a 5 exp(— a 6 p±) 


(15) 


irrespective of the angle. (Note that in the end, these details really do not matter because 
we will be developing our own parameterizations to be used in the transport codes, and 
the end results do not depend on how the Ranft parameterization is defined.) 
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Table 2. Parameters for Equations (13) and (14) or (15) and (16) 
(The numbers are the same for all reactions listed.) 


Reaction 

a x 

«2 

a 3 

a 4 

O 5 


7T + + p -4 p + X 
7 r + + p — >■ n + X 
7r~ T p — v p T X 
7r~ + p — > n + X 

0.94 

0.86 

-3.37 

3.78 

0.47 

3.6 


Substituting the Ranft parameterizations for 


drN 

dp\\dp ± 


gives 


dr a 


^"inelastic 


Eai 


[ eX P( _a 4P±) + a 5 ex p( — 06P±)] (16) 


7 o / I-. ~ inelastic ^ r~ 

dp 6 /E 2 -Ks/s 

where E = E n and p 3 , p\\ and p± are also all final state nucleon variables. 


These parameterizations are compared to data in figures 11 - 15. Clearly the Ranft 
parameterization provides a poor fit to the data. 

3. New Parameterizations of Lorentz Invariant Differential 
Cross Section for Pion-Proton Collisions 


Due to the poor fit of the Ranft parameterizations compared to the data, new pa- 
rameterizations were developed, as described below. They are compared to the Ranft 
parameterizations in figures 1 - 15. 

3.1 Pion Production 


The new parameterizations for pion production in the present work are 


and 


d 3 a 

djP/E 




2t r 


a 5 exp(— 4^p^) + a 7 


exp(-a 8 p / j P ) 
for p \ j < 0 (i.e., 9 > 90°) 


A, . ^ 

-L* ^inelastic 0 

Li T 


IE 


ai exp(— 4— p^) + a 3 exp(-a 4 p Z j P ) 


r\ / • 


(17) 


(18) 


with N = 2 and pp = 0.95. The parameters a 4 • • • a 8 are the same as listed previously in 
table 1. 
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3.2 Nucleon Production 


A major problem with the Ranft parameterization is that the value of the parameter 
<23 is negative, as shown in table 2, and sometimes leads to negative cross sections, which 
is disastrous when integrating to form spectral and total cross sections. In addition, the 
Ranft parameterization blows up for large p\ \ , which is also a disaster when integrating. 

A new parameterization using the pion parameters for 7r + production was used. It is 
the same as the new 7r + parameterization mentioned previously except that N = 2 and 
pp = 0.9. Also, where m ^ appears in the pion parameterization, m p is used here instead. 

As can be seen in figures 1-15, the new parameterizations for pion and proton 
production provide a much better fit to the data. 

4. Parameterization of Spectral Distribution 

The spectral distribution jp, is related to the Lorentz invariant differential cross section 
by references 5 and 6. 

dc T fOmax d°(T 

— = 2itp — — sin OdO (19) 

dE Jo dp 3 /E K ’ 

If is given in the cm frame, then one simply performs the numerical integration 
to obtain in the cm frame. (In the cm frame 9 ma x = 9 max = 7r -) However things 
are not this simple. The parameterizations for d '^J E of Badhwar (ref. 12), Badhwar and 
Stephens (ref. 13), and Ranft (ref. 10) are given in the cm frame, but the HZETRN code 
(refs. 1-4) requires in the lab frame. To achieve the compatibility of the frames one 
integrates over lab frame variables. One must specify E or |p| and 9 in the lab frame and 
obtain E* or |p*| and 9* in the cm frame by using equations (41) and (42). These cm 
variables are input into the formula (parameterization) for d ^ E ■ Because d ^ E is Lorentz 

invariant, the number so obtained will also be the value of in the lab frame. Thus, 
this procedure gives in the lab frame. A complication is that 9 max must be determined 
in the lab frame. This transformation is discussed in appendix D. 

4.1 Parameterizations 

The new parameterizations of the Lorentz invariant differential cross section were 
numerically integrated over angle, as described above, and then functions were fit to this 
numerical data. The following functions were found to best fit the ir production data, as 
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shown in figures 16 - 25. The resulting parameterization of the spectral distributions in 
the lab frame: 

\0.05 

— = 32.5 Log | p ! , ab ] exp [ — p*’ , 2 b ) exp[-E 0 5 ] for tt + production (20) 

/p ™ \0.05 

= 60.6 Log[pu ab ]exp[-p?^ b ] exp[-E 0 5 ] for 7r" production (21) 

Note that the cross section units are in mb/GeV and all energies are in GeV; E is the 
total energy (E=T+m) of the produced particle. Actually, the technique used was first 
to fit the total cross section, as described below, because doing so gives one a good idea 
of the p\ i a b dependence. 

The proton production spectral distributions were much more difficult to fit. The 
difficulty was in obtaining the pi i a b dependence. The best method found was to fit each 
spectral distribution separately at the pi i a b values of 2, 5, 10, 50, and 100 GeV as shown 
in figures 26 - 30. The spectral distributions for the regions in between these pi j a b values 
are well approximated by the previous pi i a b fits. For example, the spectral distribution at 
7.5 GeV is well approximated by the fitted spectral distribution at 5 GeV. Similarly, the 
spectral distribution at 75 GeV is well approximated by the fitted spectral distribution 
at 50 GeV. Thus, the parameterization developed for the proton production spectral 
distributions are 

% = exp[a2 + (b2 * E 3 ) + (c2/E) + (d2/E L5 )] for Pl lab < 2.5 GeV 

ah 

= exp[a5 + b5 * E 1 " 5 + c5/E 2 ] for p! ] ab < 7.5 GeV 

= exp[al0 + blO * E 3 + clO * (Log[E]/E)] for pi j a b < 25 GeV 

= a50 + 6 50 * E + c50 * (Log[E]/E) + (d50/E) 

+e50/£' 1 ' 5 + /50 * (Log[E]/E 2 ) + g50 * exp[— E] for p : lab < 75 GeV 

= al00 + 6100 * E + clOO * (Log[E]/E) 

+ dlOO/E + (el00/E 15 ) + fl00/E 2 for p llab > 75GeV (22) 

The parameter values are: 

a2 = -21.719014, 62 = -0.010239021, c2 = 120.40895, d2 = -110.09389, a5 = 3.8409577, 

65 = -0.75010595, c5 = -8.2109936, alO = -4.9989318, 610 = -0.00873275, 

clO = 12.995241, a50 = -0.2800469, 650 = 0.0066378301, c50 = -5.6547689, 

d50 = 27.961949, e50 = -28.032802, /50 = -9.3283952, #50 = 0.93580082, 

alOO = 0.66853886, 6100 = -0.0097222382, clOO = -16.638858, dlOO = 61.302023, 

elOO = -93.829375, and/100 ^31.86854. 

The results of the “in-between” pi lab values are shown in figures 31 - 35. It can be 
seen that the results are excellent. 


7 



5. Parameterization of Total Cross Section 


5.1 Spherical Coordinates 


The total cross section a is related to the spectral distribution cross section by 


a 


/ 


da 
dE " 


dE 



( 23 ) 


where E = \J p 2 + m 2 and p = |p|. Thus, the relation to the Lorentz invariant differential 
cross section is (ref. 6) 


a = 27t 


dd 


rPm&y 
^ Pmin 


dp 


d 3 a 


p 2 sin 9 


dp 3 / E yj'p* + m 2 


(24) 


There are two ways to calculate a. One can work in either the cm or lab frame because 
a is Lorentz invariant. 


The simplest thing is to do the entire double integration in the cm frame. In that 
case, #* iax = 7r p^ in = 0 and p^ ax is discussed in appendix C. By checking that the same 
result is obtained, one could also continue calculating in the lab frame, as was the case 
with the spectral distribution discussed previously. The value of 9 meoi j a b is discussed in 
appendix D; Pm ax , discussed in appendix C, and p* nin = 0 must both be transformed to 
the lab frame. The value of p m ; n i a b will not necessarily be zero. 

5.2 Cylindrical Coordinates 


The parameterizations of the Lorentz invariant differential cross sections of Badhwar, 
Badhwar and Stephens (refs. 12 and 13) and Ranft (ref. 10) are not only given in the 
cm frame but are also written in terms of cylindrical coordinates {p±,P\\) rather than the 
spherical coordinates (|p|,0). Thus, it is always worth trying to do the integrations in 
cylindrical rather than with spherical coordinates because the integrals may be able to be 
done analytically. Unfortunately, to obtain the spectral distribution this approach is 
impossible because E is intrinsically a spherical coordinate. A single integral over J dp\\ 
will yield -j //- , or an integral over / dp± will yield /— . but neither /— or ^ will yield || 

exactly, which is what we seek (^ can however yield ^ approximately (refs. 7 and 8)). 

Nevertheless, both spherical and cylindrical coordinate integration will yield the total 
cross section a exactly. Thus, it is an excellent check of one’s work to calculate the total 
cross section by using both coordinate systems. Obviously 
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Using equation (5), this gives 


a 


2t r 


f dp\\ j dpx 


d 3 a px 

dp 3 / E ^j p 2 + p \ + m 2 


(26) 


with E = ^ Jp ^ + p 2 ± + to 2 . 

Again, this integral can be done in either the cm or lab frames; however, given that 
the original parameterizations (refs. 10, 12, and 13) are in the cm frame, it makes sense 
to do the double integral in this frame. 


The question then arises as to the integration limits. The lower limit of p = |p| in the 
cm frame is zero, and so obviously the lower limits of p\\ and p± will also be zero. For the 
upper limit, consider the result of appendix C, where p nmx was determined in spherical 
coordinates. Given that p\\ = psinO and p± = pcosO, obviously the upper limits of p 
and p i occur when sin 9 or cos# equals unity. Thus, the upper limits of p\\ and p± will 
be the same as the spherical case, namely the exact result derived in appendix C. Thus 
the final formula, working in the cm frame, is 


fPma,x. rPma,x 

o = 2-7T J dp | J dp± 


(f'rr p 

dp 3 / E Jp 2 _|_ ^2 _|_ m 2 


(27) 


where = |p;,„| from appendix C. 


5.3 Parameterizations 


The new parameterizations of the Lorentz invariant differential cross section were 
numerically integrated over energy and angle, as described previously, and then functions 
were fit to this numerical data. The following functions were found to best fit the data, 
as shown in figures 36 - 38. 


6.5 Log[pi lab] 

for 7 r + production 

(28) 

10.1 Log[p llab ] 

for tt~ production 

(29) 

0.5 Log[pi lab]Pl'lab 

for p production 

(30) 


Note that the cross section units are in mb. 

6. Conclusions 

Ranft has provided parameterizations of Lorentz invariant differential cross sections for 
pion and nucleon production in pion-proton collisions. These parameterizations have been 
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compared to some recent data. The reactions analysed are for charged pions in the initial 
state and for both charged and neutral pions in the final state at an incident pion energy of 
40 GeV. The present work has shown some inadequacies in the Ranft parameterizations, 
and it is therefore recommended to use the parameterizations developed herein for energies 
in the 1-100 GeV range. 
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Appendix A 

Lorentz Transformations 

In the typical scenario discussed in the present paper, all variables shown subsequently 
refer to the pion; e.g., E = E^, p = p^, 6 = 9 W and so on. Variables with a star * refer 
to cm variables; while variables without a star are lab variables. Thus, for example, 9 is 
the angle of the pion in the lab frame, and 9* is the angle of the pion in the cm frame. 
Exceptions are (3 and 7 which refer to the speed of the cm frame relative to the lab. 

In this appendix, the Lorentz transformations that will be needed are summarized. 
First, define the perpendicular (transverse) and parallel momenta as 

Pj_ = p sin 9 

pu = pcos9 (31) 


where 


V = |P 


The total energy is 


E = T + m = \Jp 2 + rn 2 = Jp\ + p^ + m 2 (32) 

The Lorentz transformation of 4-momenta = (E, p) from cm to lab frame, assuming 
relative motion in the parallel direction, is given by Barnett (ref. 16): 


P . L = Pi 


E \ = ( 7 l/3\( E*\ 

P\\ ) 7 )\P\\) 

where (3 is the speed of the cm relative to the lab frame and 7 
transforms as 


(33) 


-jbp- Thus ’ p = l p 


p = ^pi + p^ = \J (p* sin 9*) 2 + ( 7 (3E* + 7 p* cos 9*) 2 (34) 

An important cjuantity is the Lorentz transformation of the pion angle. Obviously, 

p* j_ p* sin 9* 


tan 9 = — 


pjl 7 f3E* + 7 pj*j 'yPE* + 7 p* cos 9* 


(35) 
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This equation can also be written (see ref. 14) as 

t p* sin 9* 

7/34?* + 7 p* cos 6* 

E* sin 9* 

( Ei i ab + m 2 )(cos0* + a) 

with a being the ratio of the cm speed to the speed of particle 3 in the cm system 
(ref. 14), 

(3E* 
a = 

p* 

For the reaction, with a 2-body final state, 


( 36 ) 


(37) 


1 + 2 ->• 3 + 4 


there is a unique relationship between the angles and momenta of the emerging particles 
(ref. 15), so that a is also equal to (ref. 14) 

a = = / Pi iab \ s + - m\ / 3g v 

P* \Ei i a b + m 2 ) [s - (m 3 + m 4 ) 2 ][s - (ro 3 - m 4 ) 2 ] 

If three or more particles emerge in the final state, then there are no longer unique 
relationships but only limiting relationships (ref. 15). For the inclusive reaction 

1 + 2 ->• 3 + X (39) 

where X is anything, including many other particles, X = 4,5,6..., then only equation 
(37) holds. The second equality in equation (38) does not hold for the inclusive reaction 
of equation (39). (The second equality would only hold if we knew all the masses and set 
m 4 equal to the masses of all the particles X.) 


The inverse transformation from lab to cm is obtained by swapping starred and un- 
starred quantities and changing the sign of the speed (3, giving 


p transforms as 


P* = P. 



( 7 -7/3 \ / E \ 

V -70 7 ) V P\ I / 


(40) 


p* = \J (p sin 9) 2 + (—'y/3E + 7 p cos 9) 2 
The inverse angle transformation is 

tan r = = P± = psin 6 

pf, —~fj3E + 7P| —7/34? + 7 p cos 9 


(41) 

(42) 
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Appendix B 

Useful Kinematic Results 

Relation Between s and Total cm Energy 

The 4-momentum vector p ^ is 

P» = (E, P ) (43) 

with the square given by 

(p^f = E 2 - p 2 = m 2 (44) 

showing that the square of this 4- vector is invariant. For the reaction 

1 + 2 — 5- 3 + 4 

the Mandelstam s variable is defined as 

s = (Pi + P2? = (Ps + Pi? 

= (Ei + E 2) 2 — (pi + P 2) 2 = (E?y + E 4 ) 2 — (p3 + P 4) 2 (45) 

Because s is a 4- vector, then s is a relativistic invariant. The center of momentum frame 
is defined so that the total 3-momentum is zero, 

Pi + P 2 = 0 = p .3 + P 4 definition of cm frame (46) 

so that the total cm energy is 

\fs = E\ + E 2 = Es + £4 = E cm in cm frame (47) 

showing that y/s is the total energy in the cm frame. 

Relation Between s and Projectile Energy in lab Frame 

The lab (target) frame is de fin ed so that the target (particle 2) is at rest, 

P 2 = 0 definition of lab frame (48) 

Thus, in the lab frame the Mandelstam invariant s is 

s = (Ei + E 2 ) 2 - Pi 
= E 2 + 2E X E 2 + E 2 — Pi 

= ml + 2E| i ab m 2 + m 2 (49) 
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where we have used 


and 


Note that 


Ei = Ei laf, in lab frame 


/v ’2 = 'in ' 2 in lab frame 


El lab = Ti lab + nil = \j p? lab + m l 


Threshold Energy: 

Consider the reaction 


(50) 

(51) 

(52) 


1 + 2— >-3 + 4 + 5 + 6 + -- - ( 53 ) 

The Mandelstam invariant s is 

s = ( pS *+^) 2 = (K + K + ^ + ^ + -") 2 

= (Ei + E 2 f - (pt + p 2 ) 2 

= (Ez + E 4 + E 5 + E (i 4 ) 2 - (p 3 + p 4 + p 5 + p 6 4 ) 2 ( 54 ) 

The particle production threshold is defined so that all the final state particles are at rest 
in the cm frame; that is, 


P3 — P4 — P5 — P6 — • • • — 0 

Thus, at threshold, the Mandelstam invariant s is 

s = (m 3 + m 4 + m 5 + m 6 H ) 2 

Because s is a relativistic invariant, we can equate it to its lab 
(49), so that 

„ _ (m 3 + m A + m 5 + m$ H ) 2 - mj 

lab threshold “T 


(55) 

(56) 


value given in equation 
- ml 



giving the threshold kinetic energy as 


El l a b threshold 


(?n 3 + m 4 + ?7r 5 + 7776 + • • -) 2 — 777? — TOo 

— 777 1 

2777 2 

(7773 + ? 77 4 + 777.5 + 7/76 + • • -) 2 — 777 2 — ? 77 | — 2 ? 77i?772 

2t77 2 

(? 77 3 + ? 77 4 + 777.5 + 777 6 + • • -) 2 - ( 777 ! + 777 2 ) 2 

2mo 


(58) 
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or 


f I lab threshold 


(E/ m /) 2 “ ( m l + m 2) 2 

2m 2 


where E / m / denotes the sum over the masses of all final state particles. 


Speed of cm Frame 

Some other useful results for the reaction 1 + 2 — ^ anything are (ref. 16) 

Ei iab + m 2 



(59) 


(60) 


and 


0 


Pi lab 

El lab + rn 2 


7/3 


Pi lab 


(61) 

(62) 
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Appendix C 

Maximum Pion Momentum in the cm Frame 

Consider the reaction 


A + B^C + X (63) 

where A and B represent the projectile and target and X is anything. Using relativistic 
kinematics, the momentum of the C particle in the cm frame is given by (refs. 5-8 and 
15) 


D * 2 

P C 


(s + m 2 — m 2 x ) 2 — 4 srr? c 
4 s 

( s — vr? c — m \ ) 2 


as written in references (refs. 5 and 6) 


4 


°c"°x 


4s 


as written in references (refs. 7, 8, and 15) 


[s - (m c + m x ) 2 ][s - (m c - my) 2 ] 
4s 


(64) 


where Blattnig et al. (ref. 6) write sy min for m\. Also, Hughes (ref. 15) writes mx = 

Ya=i m i and E' 2 f° r s - 


In determining the maximum value of p 2 ., one demands that the X particles be pro- 
duced at rest to ensure that all available momentum goes to particle C. Also, one must 
consider the makeup of the X particles specifically. 


In previous work, pion production from proton-proton collisions was considered 
(refs. 5-8). In this reaction (with C being 7r), 

p + p — > 7r + X (65) 


X must consist of at least two nucleons due to conservation of the Baryon number 
(ref. 9). The minimum threshold reaction is thus (ref. 9) 

p + p^-TT + p + p ( 66 ) 

This results (ref. 6) in s a , min = m\ = ( m p + m P ) 2 = 4m 2 , giving (refs. 5-8) 


P 


* 2 

7r max 


(s — ?n 2 — 4?n 2 ) 2 — 16m: 


2 m 2 


4 s 

[s - ( m 7 , + 2m p ) 2 ][s - ( m ^ - 2 m p ) 2 } 
4s 


for p + p^Tr + X (67) 
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For the ir production reactions in the present work, 

7T + p — >■ 7T + X (68) 

Baryon number conservation requires at the least only one baryon to be produced in the 
final state, so that the minimum threshold reaction is (ref. 9) 

7T + p 7T + p, (69) 

which results in mx = m p . Thus, for this reaction, the maximum pion momentum is 
given by 


P 


*2 

7r max 


[s - (m w + m p f][s - (m w - m p ) 2 } 
4 s 


for 7T + p — > 7T + X (70) 


For the p production reactions in the present work, 


TT + p^p + X (71) 

Charge conservation requires X to have the same charge as the incident i r. For strong 
interactions, this must also be a 7r so that the lowest threshold reaction is (ref. 9) 

TT + p — > p + TT (72) 


Thus, mx = m w , which means that the maximum proton momentum p* p max is the same as 
the maximum pion momentum given above in (73), which is invariant under interchanging 
m, r and m p . Thus, 


*2 

P p max 


[s - {m p + 7?r 7r ) 2 ][g - (m p - ?rr 7r ) 2 ] 
4s 


for 7r + p — >p + X (73) 
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Appendix D 

Maximum Pion Angle in the lab Frame 


In the cm frame, 9* varies from 0 to i r. Hughes (ref. 15) shows that the maximum 
pion angle in the lab frame is given by 


sin m ax lab 


1 - 1/2 

\|^F/75 


(74) 


and 


V 

'in — 7 


(75) 


where V = |V| is the speed of the cm frame and v* is the speed of the pion in the cm 
frame. Therefore equation (74) becomes 


sin# 


7 T max lab 


p 


1 - V 2 

~ v 2 ~ 


T7t/^7T 

7/? 


(76) 


where we have used y 2 /3 2 = and 7 2 /3 2 = j^r- with (3 mV and /3* = v* in units where 
c = 1. Now 7^/?^ = Ip^ I /m T . Also, (ref. 16) combining j3 = F - Pla f — and 7 = Ex 

lab' m 2 V s 

gives 7/3 = Thus, 


sind 


7 r max lab 


V* Ip^I 

| Plab | ^7 r 


In carrying out integrations over the angle, one uses |p* 7rraax | for |p* 


(77) 


This result is also the same as in references 5 and 6. (However, note that |pi a b| is 
denoted as p a , and m p incorrectly appears instead of m ^ in the denominator in references 
5 and 6.) 


A much clearer discussion of the maximum angle is contained in reference 14. See 
especially the discussion on pages 402 - 403 and in figure 12.4. (The discussion occurs 
in the first edition (ref. 14) of this book but not in the two later editions (refs. 17 and 
18).) Jackson (ref. 14) does not derive the formula of equation (77), but if one plots 
figure 12.4 with specific numerical values, one can see that ^ max in figure 12.4 gives the 
same numerical answer as equation (77). Figure 12.4 is simply a plot of equations (36) 
and (37). 
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Another very clear discussion of this material is in the book by Byckling (ref. 19). In 
the cm frame, 9^ varies from 0 to 7 r. As Jackson (ref. 14) points out, just above threshold, 
a will be larger than unity, and as 6* varies from 0 to 7 r, then 9 T = 9 ^ i a b will be confined 
to some forward cone 0 < 9^ < 9 W max . As the incident energy rises above threshold, a 
will be smaller than unity and 9 W will vary all the way from 0 to 7r. 

In the reaction p + p— >-7r + p + p, analyzed previously (refs. 5 and 6) one can easily 
calculate 9 W max from equation (77) or just plot the equivalent of figure 12.4 in Jackson 
(ref. 14). One then finds that there exists a 9 7rm ax only very close to threshold. (In fact, 
it is less than | near threshold.) As soon as one gets a little bit away from threshold, 
9 tt max no longer exists and therefore 9 T = ^lab varies from 0 to ir. Given that the NASA 
radiation work considers pion production from threshold up to a very long way away from 
threshold, to all intents and purposes 9^ max is irrelevant, and one can simply consider 
that 9 n i a b varies from 0 to 7r for essentially all incident energies of interest because the 
tiny region near threshold represents only a very tiny fraction of the energies relevant to 
the space radiation problem. 

For the reactions 7r+p — > tt +p, the threshold is zero, and it turns out that 9 7Tmax never 
exists. For these reactions, ^lab varies from 0 to 7 r for all incident energies of interest. 

For proton production, the maximum angle is 

V* Ip;i 

Sill Up m ax I I 

IPlabl m P 

In carrying out integrations over the angle, one uses |p* pmax | for |p* 

For proton production, the relevant reaction is 7T + p — > p + tt, which has no threshold. 
As a result, 9 P max is also constant; however, unlike the case of pion production, this 
constant is exactly Care must be taken in numerically calculating 9 pmax when using 

equation (78). Since 9 n tnay is £ when , = 1 and 6> ax is 7 r when ■ , > 1 the 

U V / P max 2 |piabl™ P p max |pi_ab|m P ’ 

numerical round-off error can cause 9 P max to oscillate between 7 r and ^ . It is recomended 
that any code using 9 P raax put the value in explicity for the reaction 7r + p — > p + 7r rather 
than calculating it by using equation (78). 


(78) 
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Figure 8. Lorentz invariant differential cross section for 7t at 0^ m = 73°. 
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Figure 27. Fit to spectral distribution for proton production in the lab frame at p j lab = 5 GeV. (Note that the 
threshold (ref. 9) for the exclusive reaction nr + p — > p + nr is at p x ]ab tbresho i d = 0.) 



.8 
.7 
.6 
.5 

do mb 
dE ’ GeV ' 

.3 
.2 
.1 

0 5 10 15 20 

E, GeV 

Figure 28. Fit to spectral distribution for proton production in the lab frame at p l ]ab = 10 GeV. (Note that the 
threshold (ref. 9) for the exclusive reaction 7t _ + p — » p + n~ is at p x lab threshold = 0-) 


n + p — > p + X 
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Figure 29. Fit to spectral distribution for proton production in the lab frame at lab = 50 GeV. (Note that the 
threshold (ref. 9) for the exclusive reaction nr + + n~ is at p x lab tbresbo ],j = 0.) 
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Figure 30. Fit to spectral distribution for proton production in the lab frame at p 1 lab = 100 GeV. (Note that the 
threshold (ref. 9) for the exclusive reaction 7t _ + p — > p + jr - is at p b lab tbkres hoid = 0.) 


n + p — > p + X 
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Numerical integration of eqs. (17) and (18) 



Tt + p — > p + X 
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Figure 32. Parameterization of spectral distribution for proton production in the lab frame at p l j ab = 7.5 GeV. 
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Figure 33. Parameterization of spectral distribution in the lab frame at p l lab = 25 GeV. 



E, GeV 

Figure 34. Parameterization of spectral distribution for proton production in the lab frame at p l lab = 75 GeV. 
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Figure 35. Parameterization of spectral distribution in the lab frame at 1 lab = 200 GeV. 



Figure 36. Parameterization of total cross section for 7t + production (eq. (28)). 
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Figure 37. Parameterization of total cross section for k production (eq. (29] 
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